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Abstract: We describe the full characterization of the biaxial nonlinear 
crystal BiBsOg (BiBO) as a polarization entangled photon source using 
non-collinear type-II parametric down-conversion. We consider the relevant 
parameters for crystal design, such as cutting angles, polarization of the 
photons, effective nonlinearity, spatial and temporal walk-offs, crystal 
thickness, and the effect of the pump laser bandwidth. Experimental results 
showing entanglement generation with high rates and a comparison to the 
well investigated j3-BaB2C>4 (BBO) crystal are presented as well. Changing 
the down-conversion crystal of a polarization entangled photon source 
from BBO to BiBO enhances the generation rate as if the pump power was 
increased by 2.5 times. Such an improvement is currently required for the 
generation of multiphoton entangled states. 
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1. Introduction 



For more than two decades, parametric down-conversion (PDC) is a central tool for the gen- 
eration of entangled photons. This is a second order nonlinear process, where a photon from 



a pump beam splits into two photons, known as signal and idler, while conserving energy 
and momentum. The down-converted photons exhibit strong correlations in various degrees of 
freedom, such as wavelength, time of emission, polarization, momentum, and position (TJ. The 
down-conversion pair generation rate depends linearly on the pump beam power, and quadrati- 
cally on the crystal thickness and on its nonlinear coefficients |2)- Since the first demonstration 
of an efficient PDC polarization entangled photon source 0, the most commonly used nonlin- 
ear birefringent crystal for this purpose is the uniaxial crystal /3-BaB204 (BBO). The reasons 
for this are its relatively high nonlinear coefficients, high transparency, and the possibility for 
phase-matching over a broad spectral window J4][5| . 

In the last decade, many quantum optics experiments have used two consequent PDC events 
J6HH). Others have used second order events of PDC 19 HTTI . These events occur when two in- 
distinguishable pump photons split into four, during the same coherence time (or pulse duration 
for pulsed pump sources). Both approaches require high efficiency of the PDC process as suc- 
cess probability is quadratic with the single pair generation probability. Later, the third order 
PDC event as well as three consequent first order events have been used to create entangled 
states of six photons lfT2Tfl"5l . Recently, four consequent first order PDC events were used to 
demonstrate an eight photon entangled state lfl6l . 

One possibility to enhance the PDC generation probability is to use thicker nonlinear crys- 
tals. The crystal length is limited by the non-collinearity of the process, that spatially separates 
the pump beam and the down-converted photons. In addition, there is the spatial walk-off ef- 
fect between the two polarizations that degrades the entanglement quality (see Sec. 12.5b . Thus, 
higher pump intensity is required. Usually, the pump beam is generated by frequency doubling 
the radiation of a TkSapphire laser in another nonlinear crystal 1 MT61 . Reported typical in- 
tensities are above 1 W, but as the doubling crystal is damaged by the high power, it has to be 
translated continuously in order to maintain stable operation j 1 2 \ . Additionally, the pump beam 
intensity can be enhanced inside a synchronized external cavity. Such a setup has been shown 
to pump a BBO crystal with about 7 W ifTTl . 

In this work, we suggest and demonstrate the use of a novel crystal with higher nonlinear 
coefficients than BBO for the generation of polarization entangled photons. It is the mono- 
clinic biaxial BiB306 (BiBO) crystal that has been introduced [flgl and characterized lfl9l as 
a nonlinear optical crystal about a decade ago. Since then, it was used in numerous frequency 
conversion experiments (for example, see Ref. 1201 . and Refs. within). BiBO was also used 
with type-I PDC for generating photon pairs with a pulsed laser source l2D and for generating 
polarization entangled photons with a continuous pump source El . It has a very broad trans- 
parency window and its nonlinear coefficients are considerably higher than those of BBO ifTHl . 
Nevertheless, the biaxiallity introduces many differences and difficulties, compared to BBO. 

This paper is organized as follows: in Sec. [2] we present the various considerations in choos- 
ing the crystal parameters. These parameters are affected by the phase-matching angles, the 
polarization direction of the pump beam and the down-converted photons, the pump beam band- 
width, the spectral and angular properties of the down-converted photons, spatial and temporal 
walk-off effects, and the effective second order nonlinear coefficient d e ff dependence on the 
pump beam direction. Section[3]describes the experimental validation of our theoretical results 
by demonstrating and quantifying the entanglement produced by using two known configura- 
tions. 

2. Investigation of PDC parameters in BiBO 

2.1. Crystal design 

In order to lower reflections and to simplify the required calculations and alignment, it is de- 
sirable to cut the crystal facets perpendicular to the designed direction of the wave vector k/ 




Fig. 1. (a) Relative orientation of the crystallographic axes {a,} and the crystal physical 
axes {e,}. (b) Orientation of the crystal physical axes {e,} and the optical indicatrix main 
axes {e® }. 4> is the angle of orientational dispersion of the principal axes, (c) The propaga- 
tion direction of the pump beam T inside the crystal is defined using the two angles iff and 
p in the wavelength independent {e,} system. 



of the fundamental (pump) wave. The phase-matching calculation for the k/ direction scans 
a quadrant of space, according to the monoclinic symmetry of BiBO. In order to choose the 
optimal phase-matching direction, the effective second order nonlinear coefficient d e ff is cal- 
culated for each direction. As an approximation for d e tf, we use the effective nonlinear 
coefficient of collinear second-harmonic generation, d^f , calculated for any direction, even 
though the phase-matching condition is not fulfilled. The optimum direction of kf within the 
range of the highest values of d e ff should allow the two cones to intersect at 90°, which is 
optimal for the photon collection efficiencies. For the selected direction of k^ as well as for 
the down-converted photons at the intersection points of the emission cones, the polarization 
orientation is calculated. Finally, the temporal and the spatial walk-offs are calculated for the 
chosen crystal parameters. 

2.2. Phase-matching calculation 

In order to find the spatial distribution of the cones of down-converted photons, we calculated 
numerically the non-collinear type-II PDC process in BiBO. We are interested in the degen- 
erated case in which the down-converted photons share the same wavelength. The most basic 
reference system that we use is the crystal physical Cartesian system {e,}. It is linked to the 
crystallographic system {a,} (see Ref. l23ll ) by \\ a^,e2 || «2 II 2-fold axis, e\ = x e$, see 
Fig. 0I a )- The point group symmetry 2 of the monoclinic BiBO crystal structure allows the 
occurrence of enantiomorphic (i.e., "left-handed" and "right-handed") species. All our samples 
for optical investigations were prepared using crystals that were grown as descendants from the 
same parent crystal and therefore posses the same handedness. For our crystals, the positive 
direction of «2 ( an d e i) corresponds to a positive sign of the pyroelectric coefficient (at 
constant stress) and to a negative sign of the longitudinal piezoelectric coefficient if 222 112411251 . 

In BiBO, the principal axes {e?} of the optical indicatrix coincide with the {e,} system 
only for e\ = ei while e\ and e\ change their orientation with wavelength. This orientational 
dispersion is illustrated by the angle <t> = L (63,63) in Fig.OJb). For the fundamental and the 
down-converted wavelengths used in this work (A/ = 390 nm and X c j c = 780 nm), 4> equals 
43.8° and 46.9°, respectively fl9l . 

Our calculations of the collinear and non-collinear PDC phase-matching cases 12611 . basically 
follows the calculation strategy described by Ref. l27l . For a chosen direction T of the funda- 
mental wave vector k/, we define the propagation direction in spherical coordinates ( y/, p ) with 
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Fig. 2. (a) Stereographic projection of collinear type-II phase-matching angles for 
Ay =390nm (black line), with several non-collinear down-converted circles for different 
propagation directions of the fundamental wave. The chosen working point for this work 
is marked with X. (b) Stereographic projection of PDC in BiBO for 1/a(T) = 63.5° and 
p(T) = 53.5°. The vectors P and R are indicated, as defined in the text. 



respect to {e,} (see Fig.QJc)). The phase-matching conditions are satisfied when 

Ak = k sigm i +k mer -k f = 0, ( 1 ) 

where Ak is the phase-mismatch vector, and k s i gna i and kuii er are the wave vectors of the down- 
converted waves. First, we find the collinear phase-matching angles, as in this case Eq. Q]be- 
comes scalar and simple to solve. Then, we use a search algorithm around the collinear direc- 
tion to find the non-collinear directions that correspond to the minimal values of Ak. We have 
chosen a numerical threshold value of |^| < 5 x 10~ 5 . Photons are emitted into two cones 

with different, and not necessary perpendicular, polarizations. The stereographic projections of 
several down-converted emission cones onto the (e\,ej,) plane are presented in Fig.[2|a). This 
projection preserves angles and projects circles in three dimensions as circles on the plane (26| . 
Each two tangent circles represent a non-collinear solution, where the direction of the funda- 
mental wave k/ is their collinear intersection point. The down-converted photons experience 
refraction when they emerge from the crystal to air, which depends on their propagation di- 
rection and their polarization. The calculation results given in this work are of the photon's 
properties outside the crystal. For our wavelength parameters, the phase-matching calculations 
resulted in a suitable direction T with spherical coordinates y/ = 63.5° and p = 53.5°. In this 
case, the two down-converted cones intersect at an angle of 90° and the intersection points 
are separated by 6.9 ± 0.2°. In order to simplify the crystal alignment process, it is convenient 
to define a sample reference system according to the PDC emission results. The direction of 
the wave vector k/ of the fundamental wave is parallel to T. T is also normal to the input 
facet of the sample. We define P to be the vector connecting the two cones intersection points 
(see Fig. |2b)) and R the vector that connects the most distant points on each circle. Conse- 
quently, T, P, and R form an orthogonal set. The BiBO samples used in our PDC experiments 
have spherical coordinates (>,p) of T = (63.5°, 53.5°), P= (-80.6 ± 0.1°, +30.9 ± 0.1°), 
and R = (— 1.4±0.1°, — 17.4±0.1°). The errors result from the finite grid resolution of the 
calculation for the intersection points. 
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Fig. 3. (a) The relevant photon polarization directions for the designed crystal. The elliptical 
cross-sections of the wavelength dependant indicatrix are marked for the fundamental beam 
and the two cones intersection directions. Note that the ellipticity of the cross-sections is 
exaggerated for clarity reasons. The long and short semi-axes of the cross-sections indicate 
the polarization directions of the slow and fast waves, respectively, (b) Experimental results 
of the normalized intensity of the down-converted photons as a function of the fundamental 
wave polarization angle. 



2. 3. The photons ' polarization 

For any light propagation direction inside a non-cubic crystal, there are two orthogonal modes 
of the dielectric displacement field Di and D2, each with a different corresponding refractive 
index. In uniaxial crystals, such as BBO, these two modes are known as the ordinary (o) and 
extraordinary (e) polarizations, while in biaxial crystals they are known as the fast (f) and slow 
(s) polarizations, both behaving in general as an extraordinary wave |28l - 

When choosing the crystal parameters, we need to consider the polarization of the pump 
beam and the down-converted photons. It is possible to calculate the directions of the dielectric 
displacement vectors Di and D2 of the two linearly polarized waves in respect to the physical 
axes {e,}. However, it is more convenient to define the photon polarizations with respect to 
the P and R directions. In the BBO crystal, the pump beam is polarized along the R direction, 
one cone is polarized in the same direction, and the other cone is polarized in the P direction. 
In BiBO, the pump beam should be polarized in its fast polarization mode in order to achieve 
maximal conversion efficiency, which usually differs from these convenient directions. For the 
general case, we define the cartesian coordinates of the propagation direction T by the unit 
vector (x,y,z) in the optical indicatrix system {e?}. Using the Sellmeier formula for BiBO |fl9l , 
we calculated the wavelength dependant principal refractive indices [n x < n y < n z ). From them 
we derived the slow and fast refractive indices (27). Using these refractive indices, the ratios 
between the components of the normal polarization modes (i.e., the components of the unit 
vectors along the displacement field vectors D,-) are given by l28l 

n 2 x x n\y n\z 

(nf-ni) (nf-rf) (nf-n 2 z ) 

where i stands for 'fast' or 'slow'. For our crystal parameters, the fast polarization mode of the 
pump beam was calculated to be 13. 2 ±0.1° from P, as shown in Fig. [31 a). We also measured 
this value by rotating the pump polarization direction with a half-wave plate. At each rotation 




Fig. 4. (a) The calculated d e ft [pm/V] of BiBO. Each contour line marks a step of 
0. 14pm/V. The thick black line represents collinear type-II phase-matching directions for 
Xj =390 nm. The X symbol marks the chosen pump direction in this work, (b) The Spatial 
walk-off angle [deg.] for BiBO. Each contour line marks a step of 0.45°. 



step we took a picture of the down-converted circles. In Fig. [3fb) we plot the normalized in- 
tensity of the down-converted photons vs the polarization angle. Setting 0° parallel to P, the 
maximal value was obtained at an angle of 11.6 ±0.3° from P, within the crystal fabrication 
errors. 

In order to calculate the polarization of the down-converted photons at the cones intersection 
points, we need to consider these two propagation direction inside the crystal. For the down- 
converted photons propagating at the top left (bottom right) intersection point in Fig. 0a), the 
polarization direction of the fast wave is at 14.5° (16.0°) from P. The slow polarization modes 
are perpendicular to the fast modes. 

2.4. The effective second order nonlinearity 

The effective strength of the second order nonlinear coefficient d e ff is an important consid- 
eration for the crystal design. For the uniaxial BBO crystal, there is an analytical expression 
that appears in Ref. (29). Using the BBO d matrix elements from Ref. Q and our wavelength 
parameters, a maximal value of d e ft = 1.1 5 pm/V is calculated. The calculation assumes a 
collinear type-II phase-matching process. 

A rigorous treatment of biaxial crystals appears in Ref. |30l . We used the relevant formula 
for dgjj of collinear type-II phase-matching in the {e®} reference system (the fsf indices re- 
fer to the pump and the down-converted photon polarization modes). The calculation results 
were rotated to the {e,} reference system, where for the parameters used in this work we get 
d e ff = 2.00pm/V. The calculation considers a wavelength of 780 nm, although there is almost 
no wavelength dependency. For this calculation we used the four d matrix elements given in 
Ref. Il30ll . The results for any T direction are shown in Fig.Ua). Note that because the calcula- 
tion assumes collinear propagation, the results have significant meaning mainly in the vicinity 
of the collinear phase-matching curve. Furthermore, we have also removed the Kleinman sym- 
metry assumption of Ref. 13011 and derived a formula containing the eight d matrix elements 
given in Ref. BTI . This generalization resulted with a similar value (d e ff = 2.02 pm/V). The 
almost doubled value of the nonlinear parameter of BiBO compared to BBO promises a major 
advantage for the generation of entangled photons. 



2.5. The spatial walk-off angle 

During the propagation through a birefringent crystal, the Poynting vector may point away 
from the direction defined by the k vector, depending on the beam polarization [32]. This phe- 
nomenon is called spatial walk-off. It should be taken into consideration when designing a 
polarization entangled photon source since it can create spatial labeling of the down-converted 
photons, which in turn will reduce the entanglement quality. The spatial walk-off angle SWO 
between the Poynting vector and the k vector, together with the crystal thickness L, determines 
the overall spatial walk-off. The pump beam spot-size at the crystal should be large compared 
to the spatial walk-off in order to prevent the labeling effect J3). 

In uniaxial crystals, such as BBO, an ordinary photon's k vector and Poynting vector have 
the same direction while an extraordinary polarized photon deviates from that direction by an 
angle that can be calculated using a simple analytical expression l33l . In biaxial crystals, such 
as BiBO, both the fast and slow polarized photons deviate from the direction defined by the 
k vector while passing through the crystal. The spatial walk-off angle in this case is the angle 
between the two down-converted photons' Poynting vectors. 

We present here the results of a numerical approach for the walk-off calculation for BiBO. 
The direction of the Poynting vectors of the slow and fast down-converted photons are normal to 
the surface of the corresponding indicatrix. For each photon, we calculated three wave vectors 
with small deviations from their propagation direction k. We then found the plane that contains 
these three vectors. The direction normal to this plane is the direction of the Poynting vector. 
The angle between the two Poynting vectors of the slow and fast photons is the required walk- 
off angle. Note that it is also possible to treat this problem analytically, but as our numerical 
results are sufficiently accurate, we leave the rigorous treatment for a later work. 

We calculated numerically the spatial walk-off angle in BBO and BiBO for a wavelength 
of Xdc = 780 nm. We have validated our numerical approach by comparing its results to the 
analytical expression for BBO 0331 . The typical deviation between the numerical and analytical 
calculations is about 10~ 6 degree. For collinear PDC in BBO the walk-off angle is 9 swo = 4.15°, 
corresponding to an overall walk-off of 145 fim for a 2 mm thick crystal. The results for BiBO 
are presented in Fig.Ub). For our crystal parameters, the calculated walk-off values are Q swo = 
3.55° for one of the cones' intersection points and 9 swo = 3.6° for the other. These results 
correspond to a deviation of about 95 jJ.m for the 1.5 mm thick crystal used in our experiments. 

2.6. The temporal walk-off 

As their name suggests, the two polarization modes propagate through the birefringent crystal 
with different group velocities. This may cause temporal distinguishability between the slow 
and fast photons. This phenomena is known as temporal walk-off. A birefringent crystal of 
thickness L separates the photons by 

L L ,nl n( , An r 

8T = =£(-£-— )=L—, (3) 

v s Vf c c c 

where c is the speed of light in vacuum, and v s (v/) and n s r (n() are the group velocity and the 
ray refractive index of the slow (fast) photon, respectively. The ray refractive index «, and the 
refractive index n are related via n r = ncos a, where a is the angle between the k vector and the 
corresponding Poynting vector (32). The problem is more significant when 8T is comparable 
to or larger than the coherence time z c . We addressed this issue with two methods. The first 
is to add two compensating crystals, cut at the same directions as the generating crystal but of 
half the thickness, in each down-conversion path Q. The second approach is to overlap the two 
photons at a polarizing beam splitter (PBS) 1341 , as will be described later in Sec. 13. II 




Fig. 5. (a) Down-converted photon circles through a 3nm bandpass filter, from a 2 mm 
thick BBO crystal, (b) Down-converted photon circles through a 3 nm bandpass filter, from 
a 2.7 mm thick BiBO crystal with different pump wavelengths, as indicated. Several lower 
circles are cropped due to the filter size. 
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Fig. 6. (a) Stereographic projection of non-collinear type-II PDC processes in BiBO 
with different pump wavelengths. The two inner circles originate from a pump beam of 
A = 389 nm, while the two outer circles from A = 391 nm. The thicker ring (left, red) is 
polarized slow, while the thinner one (right, green) is polarized fast. The X symbol marks 
the pump direction for the BiBO crystal in this work, (b) Stereographic projection of non- 
collinear type-II PDC processes in BiBO with a pump wavelength of 390 nm and different 
down-converted wavelengths. The two inner circles wavelength is 781.51 nm (left, red) and 
778.5 nm (right, green) and the two outer circles are of the opposite process. 



For our crystal parameters, An r is approximately 0.05 for BBO and 0.15 for BiBO, which 
results with 8T = 330 fs for a 2 mm thick BBO and 8T = 750 fs for a 1.5 mm thick BiBO. 
Compensation is required in both cases as these values are larger than x c = 180 fs, the coherence 
time that corresponds to the used 3 nm filters. 

2. 7. Pump bandwidth and the entanglement quality 

One advantage of down-converting a pulsed source over a continuous source is its energy con- 
centration in a short coherence length which increases the probability of higher order PDC 
events. Furthermore, its timing information is inherited by the down-converted photons. How- 
ever, the pulses broadband spectrum can cause a variety of undesired effects that decrease the 
entanglement quality. 
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Fig. 7. A comparison of the measured (open circles) and calculated (solid circles) normal- 
ized down-converted circles radii of the slow (dashed red) and fast (solid blue) photons. 



Table 1. |£ for X = 780 nm in nm" 1 



BBO BiBO 
slow 3.15 x 1(T 5 7.0 x 1(T 5 
fast 2.85 x 1(T 5 5.0 x 1(T 5 



Figure [3a) presents a picture of the down-converted photons from BBO recorded by a sen- 
sitive CCD camera through a 3 nm bandpass filter. The pump wavelength is A/ = 390 nm with 
a full width at half-maximum (FWHM) of ^ 2 nm. Figure |5jb) shows pictures of the down- 
converted photons from BiBO for a few pump wavelengths. From these pictures we note a 
clear difference between the widths of the two BiBO circles, which is much smaller for BBO. 
For BiBO, the lower circle, which is made out of the slow polarized photons, changes much 
more than the fast polarized circle. This is due to the difference in the dispersion of the refrac- 
tive indices at a wavelength of 780 nm for slow and fast polarized photons in this propagation 
direction. The calculated dispersion values from the Sellmeier formulas for BBO and BiBO are 
presented in Table Q] In the BiBO case there is a 40% difference, while in BBO the dispersions 
differs only by about 10%. To ascertain these results we calculated the processes that corre- 
spond to down-conversion of wavelengths at the FWHM values of the pump 2 nm spectrum 

fast(389nm) — > slow(780nm) + fast (776.01 nm), 
fast(389nm) — ► s/ow(776.01rcm) +fast(780nm), 
fast (391 nm) — > slow(780nm) + fast (784.01 nm), 
fast(39l nm) — ► s/ow(784.01 nm) + fast (780 nm). 

We present on a stereographic projection only the circles of X c / C = 780 nm (Fig. |6ja)). The 
circles angular radii are measured and normalized by the radius of 780 nm circles from down- 
converting 390 nm photons. Figure|7]presents a comparison between the numerically calculated 
radii and those measured from Fig. |5jb). The calculated (measured) slopes for the two polar- 
izations differ by a factor of 3.65 ± 0.15 (2.55 ± 0.05). The calculated slow polarization circle 
is thicker than the fast polarization circle by 2.8 ± 0. 1 times. 

In order to separate the effect of the pump bandwidth from the effect of the filter bandwidth, 




Fig. 8. The BiBO (a) and BBO (b) collinear type-II PDC spectra. For BiBO (BBO), the 
spectrum of the fast (extraordinary) photons is presented by a solid blue line, while that 
of the slow (ordinary) photons' by a dashed red line. In both cases, the crystals' thickness 
is 2 mm, the filter bandwidth is 3 nm, and the pump bandwidth is 2 nm. Spectral overlap 
is 89.6% for BiBO and 98.2% for BBO. Insets: Phase-matching spectral dependency be- 
tween the slow (ordinary) and the fast (extraordinary) photons from BiBO (BBO). The 
spectra aspect ratios are 1:3 and 2:3 for BiBO and BBO, respectively. 
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Fig. 9. (a) Spectra overlap as a function of the crystal thickness with a 3 nm bandpass filter 
for BiBO (blue squares, solid line) and BBO (red circles, dashed line) crystals, (b) Spectra 
overlap as a function of the filter bandwidth for 2 mm thick BiBO (blue squares, solid line) 
and BBO (red circles, dashed line) crystals. 



we have also calculated the circle widths due to the filter bandwidth for a 390 nm pump. We 
consider the processes that result with photons at the FWHM of the 3 nm filters 

fast (390 nm) — > slo w(778 .5 nm) + fast (78 1 .5 1 nm) , 
fast(390nm) — ► s/ow(781.51 nm) + fast (778.5 nm). 

The results are presented in Fig. 1Mb). There is no significant effect due to the filter's width. 
Thus, the slow polarized circle larger width is attributed to its higher dispersion, that results 
with the asymmetry shown in Fig. [2b). The filters bandwidth do not add asymmetry between 
the circles. This conclusion suggests that a symmetric PDC picture may be obtained from BiBO 
using a continuous pump source. 
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Fig. 10. The experimental setup. See text for details. 



In order to evaluate the effect of these results on the quality of the generated entangled state, 
we calculated the spectra of the down-converted photons form a pulsed source in BiBO and 
in BBO. Our calculation was based on the work of Grice et al, that was previously applied to 
BBO 1351 . As before, we used the collinear approximation. The normalized overlap between 
the two down-converted photons' spectra corresponds to the quantum state visibility. The pump 
bandwidth and the crystal thickness also influence the down-converted spectra and thus, should 
be considered when designing such a polarization entangled photon source. Bandpass filters 
with the proper bandwidth can enhance the overlap between the two down-converted photons, 
and thus reduce the distinguishability between them. Figure[8]presents calculations of the down- 
converted spectra for a 2 mm thick BiBO and BBO crystals, assuming spectra with a FWHM 
of 2 nm for the pump photons and with 3 nm for the bandpass filters. The overlap between the 
integrated spectra of the two photons for BiBO and BBO are 89.6% and 98.2%, respectively. 
For the entanglement measurements, we used a 1.5 mm thick BiBO crystal, with a calculated 
spectral overlap of 92.8%. Figure |9ja) presents the dependency of the BiBO and BBO spectra 
overlap on the crystal thickness for a 3 nm filter. The same spectra overlap as a function of the 
filter bandwidth for a 2 mm thick crystal is presented in Fig. 0b). Although it seems as BBO 
can perform better than BiBO, spectral distinguishability can be eliminated 11341 . 

3. Entanglement measurements 

3.1. The experimental setup 

The setup used in this experiment is presented in Fig. |T0j The radiation of a mode-locked 
Tksapphire laser at 780 nm is up-converted to 390 nm by second-harmonic generation (SHG). 
The beam is focused by a lens (L) on the BBO or BiBO crystal (NL). The spatial modes of 
the pump beam and the down-converted photons are matched to optimize the collection effi- 
ciency 11361 . The photons are coupled into single mode fibers (SM), where their polarization 
is adjusted by polarization controllers. The relative propagation delay between the two optical 
paths is adjusted by translating one of the fiber ends with a linear motor (M). A quarter-wave 
plate (QWP) and a half-wave plate at each path are used for the quantum state tomography. 
The photons are spectrally filtered by using 3 nm wide bandpass filters (F) and coupled into 
multimode fibers that guide them to the single-photon detectors (SPD). 

We tried two configurations in order to remove the temporal and spectral distinguishability of 
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Fig. 1 1 . Results with configuration I. (a) Visibilities vs the twofold coincidence rate in three 
polarization bases: HV (black squares), PM (red circles), and RL (blue triangles). Straight 
lines represent linear fits, calculated without the last three points, where stimulation is more 
significant, (b) Twofold coincidence rates vs pump power. The solid black line represents 
the quadratic fit and the dashed red line the linear slope at low pump powers. 
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Fig. 12. Results with configuration II. (a) Visibilities vs the twofold coincidence rate in 
three polarization bases: HV (black squares), PM (red circles), and RL (blue triangles). 
Straight lines represent linear fits, (b) Twofold coincidence rates as a function of the optical 
path difference. The red circles correspond to a projection to the \<j> + ) state and the black 
squares, a projection to the \ <p~) state. Blue triangles represent coincidence events from the 
same side. 



the down-converted photons. The elements used in each configuration are labeled I and II in Fig. 
[TOl In the first configuration (I), the photon polarizations are 90° rotated by a half-wave plate 
(HWP), and temporal and spatial walk-offs are corrected by compensating crystals (CC) of half 
the thickness of the generating crystal. In the second configuration (II), two perpendicularly 
oriented Calcite crystals (arrows indicate the optical axis direction) are used for aligning the 
birefringent phase (BP). The photons are then overlapped at a PBS. 

3.2. Experimental results 

We generated polarization entangled states with a 1.5 mm thick BiBO crystal and compensated 
for distinguishability effects with two configurations 1131 1341 (see Sec. 13. 11 1. In order to character- 
ize the entanglement quality, we recorded visibilities J5) at three polarization bases (horizontal 
and vertical linear polarizations (HV), plus and minus 45° linear polarizations (PM), and right 
and left circular polarizations (RL)). Full quantum state tomography was also performed. Com- 
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Fig. 13. Real parts of the measured density matrices for the two configurations. Imaginary 
values are smaller than 0.08 and therefore not presented, (a) Configuration 1, 40 mW pump, 
(b) Configuration I, 300 mW pump, (c) Configuration II, 42 mW pump, (d) Configuration 
II, 3 10 mW pump. 

parison is made with results obtained using a 2 mm thick BBO crystal, in a setup optimized for 
its parameters. 

Using configuration 1, we generated the Bell state. The recorded visibilities were 

V H v = 93.4 ±0.8%, V PM = 87.5 ±0.6%, and V RL = 85.8 ±0.5%, with typical twofold co- 
incidence rates of 34000 Hz (pump power of 300 mW). The detection efficiency, i.e., the ra- 
tio of the rates of coincidence events and of single events, was 8.7 ±0.7%. By lowering the 
pump power with a variable neutral density filter these values were improved (see Fig. [TO a)). 
The visibilities dependence on the pump power is due to higher order PDC events that result 
from the improved generation efficiency. The extrapolated visibilities at zero pump power are 
V HV = 98.5 ± 0.1%, V PM = 92 ±0.1%, and V RL = 89.9 ±0.1%. Another manifestation of the 
high generation efficiency is the stimulation process that can be seen from the dependence of 
the twofold coincidence rate on the pump power (see Fig.fTTTb)). A quadratic function fits the 
data well, and clearly deviates from the linear slope at low pump powers, a clear signature of 
stimulated PDC. Density matrices were measured both at low power (P=40mW, Fig. [T3l a)) 
and high power (P=300mW, Fig.[Hb)) El. Their fidelities are 0.91 ±0.01 and 0.88±0.01, 
respectively, calculated using a maximal likelihood algorithm 11381 . 

The second compensation scheme we have used (configuration II), was first suggested and 
experimentally demonstrated by Kim et al. f34l . Using this configuration, we generated the 
|0+) Bell state. The recorded visibilities were V H v = 90.2 ± 0.8%, V PM = 86.9 ±0.9%, and 
Vrl = 85.9 ±0.9%, with typical twofold coincidence rates of 22500 Hz (pump power of 
310mW, see Fig. fT27 a)). The lower count rates can be attributed to the lack of spatial walk- 
off compensation due to the missing compensating crystals. The detection efficiency was 
7.2±0.5%. The extrapolated visibilities at zero pump power are Vhv = 96.9 ±0.1%, Vpm = 
93.3 ±0.1%, and Vrl = 90.9 ±0.1%. The density matrices for pump powers of 42 mW and 
310mW are shown in Figs, \l3\c) and [T3ld), respectively. Their corresponding calculated fi- 
delities are 0.94 ±0.01 and 0.90 ±0.01. 

There is another way to evaluate the entanglement quality when the second configuration 
is used. Scanning the path difference before the PBS and recording coincidences at the 45° 
rotated base, simultaneously projects on the |0 + ) and states (see Fig. [T27 b). pump power 
is 320 mW). The dip visibility is Vpm = 86.4 ±0.5%, similar to the high pump power visibility at 
the PM basis. This value is affected by contributions from high order events. The contribution 
of the second order term to the visibility can be estimated from the coincidence rate of two 
orthogonally polarized photons at the same PBS output port. After subtracting the second order 
events contribution, the visibility becomes Vpm = 91.1 ±0.6%, which is comparable to the PM 
visibility at low pump power. 

For comparison, we have generated polarization entangled photons in the \<p + ) Bell state 



from a 2 mm thick BBO crystal with configuration I. The measured visibilities in the three 
polarization bases Vhv,Vpm, and Vrl were 95 ± 1%, 91 ± 1%, and 89.5 ±1%, respectively. The 
typical twofold coincidence rate was 37500 Hz for a pump power of 410 mW and detection 
efficiency of 13 ± 1%. We have also measured the density matrix and calculated the state fidelity 
to be 0.95 ±0.01. 

A meaningful comparison between BiBO and BBO should take into account the differences 
between the two crystals we checked. The two crystals were also measured in different setups, 
but these were individually optimized to optimize the collection efficiency, which depends on 
the crystal parameters. The BBO crystal we used was 2 mm thick and anti-reflection coated, 
while the BiBO was only 1.5 mm thick and uncoated. The thickness difference accounts for a 
factor of 1 .78, as the PDC efficiency depends quadratically on the crystal thickness J2) . The lack 
of coating for BiBO also accounts for a 12 ± 1% loss, assuming the pump beam and the down- 
converted photons hit the crystal facets perpendicularly. It should also be considered that, due 
to some technical issues, we pumped the two crystals with different powers. Thus, we calculate 
the down-conversion efficiency as the number of detected pairs per second, per mW of pump 
power, per mm 2 of crystal thickness. The efficiency values for BiBO and BBO, as measured 
in configuration I, are 58 ± 1 Hz mW^'mm -2 and 23 ± 1 Hz mW _1 mm~ 2 , respectively. These 
values account for an improvement by 2.5 ± 0.15, compared to the 3.09 ratio predicted by the 
calculated d e ff values of BiBO and BBO (see Sec. 12.41 ). 

4. Conclusions 

We have studied the various properties of the biaxial BiBO crystal, which are relevant for 
utilizing it as a polarization entangled photon source using non-collinear type-II PDC and a 
pulsed pump source. Theoretical and numerical treatment of the relevant crystal parameters is 
presented. We calculated the crystal cutting angles, the polarization directions, temporal and 
spatial walk-offs, and the effective nonlinear coefficient. We have also demonstrated the effects 
of crystal dispersion and the broad spectrum of the pulsed pump on the angular and spectral 
properties of the down-converted photons, and therefore on the entanglement quality. The ex- 
perimental results demonstrate the higher efficiency of BiBO compared to the commonly used 
BBO, and the potential BiBO has as an ultra bright source of entangled photons. Although it 
focuses on BiBO, our work can be considered as general guidelines for considering any other 
biaxial nonlinear crystal as a non-collinear type-II polarization entangled photon source. As 
there are a growing number of quantum optics experiments that require highly efficient PDC 
sources, we hope that this work will encourage the use of BiBO as a source for polarization 
entangled photons. 
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